Cascading Gravity is Ghost Free 
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Abstract 

We perform a full perturbative stability analysis of the 6D cascading gravity model 
in the presence of 3-brane tension. We demonstrate that for sufficiently large tension 



on the (flat) 3-brane, there are no ghosts at the perturbative level, consistent with 



results that had previously only been obtained in a specific 5D decoupling limit. 



C*~> ■ These results establish the cascading gravity framework as a consistent infrared 

> ; 

. modification of gravity. 

o 

i-H 

(N ■ I. INTRODUCTION 

o 
o 

Ten years after its discovery, cosmic acceleration remains one of the central puzzles in 
cosmology. Although the data is thus far consistent with dark energy plus General Relativity 
(GR), a tantalizing alternative is that new gravitational degrees of freedom on cosmological 
scales are responsible for late-time acceleration. The most compelling idea along these lines is 
that the graviton has a small mass or width, of the order of today's Hubble parameter, which 
could account for the apparent smallness of the cosmological constant via the degravitation 
mechanism [1-3]. 

The central difficulty in constructing any consistent infrared (IR) modified theory of 
gravity is the avoidance of ghosts. For instance, massive gravity theories with a hard mass 
for the graviton [4] are well-known to be unstable [5-7]. More promising constructions 
rely on branes and extra dimensions [8-15], such as the Dvali-Gabadadze-Porrati (DGP) 
model [8]. The normal branch of the DGP model is perturbatively ghost free, in contrast to 



2 



the self-accelerating branch [16, 17], and thus represents a perturbatively consistent infrared 
modification of gravity in which the gravity has a soft mass, i.e., it is a resonance. 

The Cascading Gravity framework [9-13] proposed recently generalizes the DGP model 
to higher dimensions. There are two motivations for considering this broader set-up. As 
theoretical motivation, it offers a new class of IR modified gravity theories whose properties 
are therefore worth exploring. In particular, the lower momenta dependance of the mass 
could lead to important consequences for degravitation where the fundamental cosmological 
constant, here arising with tension of the 3-brane, could be large but give rise only to a 
small backreaction on the geometry on timescales of the order of the graviton Compton 
wavelength. As observational motivation, the model carries distinguishable signatures that 
could be observable at late times [18]. 

In Cascading Gravity, our 3-brane is embedded in a succession of higher-dimensional 
branes, each with their own intrinsic Einstein-Hilbert term. In the simplest realization, with 
6D bulk space-time, the action is 

M i f 6 Mf f 5 

Sqd cascading = \ d xy/—g 6 R 6 + — - / d x^—g^R^ 

./bulk ^ </4-brane 

+ / d^v^ ( ^$R 4 + 40 . (1) 

J 3-brane V Z J 

There are two characteristic crossover scales, m 5 = M|/M| and m 6 = Mg/M|. Assuming 
rriQ <C m 5 , the gravitational potential on our brane cascades from a 1/r (4D gravity) regime 
at short distances, to a 1/r 2 (5D gravity) regime at intermediate distances, and finally 
to a 1/r 3 (6D gravity) regime at large distances. Remarkably, the codimension-1 kinetic 
term makes the 4D propagator finite [9, 10], thereby regulating the logarithmic divergence 
characteristic of pure codimension-2 branes [19, 20]. Some of the cosmological implications 
of this model have been explored in [18, 21, 22]. 

The ghost issue is trickier. Perturbing around the flat space solution with empty branes, 
one finds that a ghost scalar mode propagates. This is seen most directly from the tensor 
structure of the one-graviton exchange amplitude between conserved sources in the UV limit: 

•4 ~ Tr ■ ■ - \n^v) - I T, ■ ^- ■ Tl . (2) 

We have conveniently separated terms as the sum of a massive spin-2 contribution, with the 
well-known 1/3 coefficient, plus a contribution from a conformally-coupled scalar tc. The 
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problem is with this scalar mode. The last term in (2) is negative, indicating that n is a 
ghost. This UV behavior is identical to other higher- dimensional scenarios [23, 24]. 

However, it was immediately noticed [9] that the ghost is removed if the codimension-2 
is endowed with sufficiently large tension A. (In analogy with a cosmic string in ordinary 
4D gravity, adding tension to a codimension-2 defect leaves the induced geometry flat but 
creates a deficit angle in the extra dimensions.) Instead of (2) the exchange amplitude now 
depends on the tension A: 

A ~ ^•^•( 7 ^-^) + rfir — 7\ % -¥hs- r: - (3) 

The 7r contribution is healthy, and the ghost is absent, provided that 

A > \m\Ml . (4) 

Note the essential role of the codimension-1 brane. In the limit M 5 — > where its induced 
gravity term disappears, the lower bound in (4) diverges. In particular, adding tension 
cannot cure the ghost in the pure codimension-2 DGP case [23, 24]. 

A limitation of the calculation presented in [9] is that (3) and (4) were obtained through 
a decoupling limit of (1), namely M 5 , M 6 — » oo keeping A 6 = (mgMf ^Z 7 fixed. (The deriva- 
tion of these results using the decoupling limit is reviewed in Appendix A.) Although the 
decoupling limit in DGP has been shown to faithfully reproduce much of the phenomenology 
of the full-fledged higher-dimensional theory [16, 25], a complete 6D calculation is clearly 
needed to establish unequivocally the consistency of the cascading gravity framework. In 
particular, the analysis of [9] could not demonstrate the absence of ghost modes that disap- 
pear in the decoupling limit. 

In this paper, we study the issue of stability rigorously by perturbing the full 6D action (1) 
around a background including tension on the 3-brane. The background geometry is flat 
everywhere, but the extra dimensions show a deficit angle due to the codimension-2 source, 
with the 3-brane located at the conical singularity. 

We are immediately faced with a puzzle: brane tension, by virtue of being a source for 
gravity, cannot by itself modify the graviton kinetic term. Indeed, the 3-brane part of the 
quadratic lagrangian density takes the initial form 

^3-branc = J d'x 7rD 4 7T - 2^) , (5) 
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where tt is the scalar perturbation of the 4D metric. Thus it acquires a localized mass term 
from A, but as expected its kinetic term is unscathed and remains negative. Incidentally, 
even the mass term is puzzling — n is a Goldstone mode and should enjoy a shift symmetry. 
How can (5) be consistent with the conclusions of the decoupling analysis? 

The resolution lies in the other scalar modes: integrating out the non-dynamical degrees 
of freedom in the bulk results in contributions localized on the 3-brane, which generate 
A-dependent additions to the 7r kinetic term. The final 3-brane lagrangian for n is 



Provided that A > 2m 2 1 M|/3, the kinetic term for n is positive, in precise agreement with the 
decoupling limit. Meanwhile, the puzzling mass term gets canceled by a bulk contribution. 
Our analysis shows unambiguously that the cascading framework is ghost-free, at least 
perturbatively, provided that the codimension-2 brane is endowed with sufficiently large 
tension. Whether the theory is stable non-linearly remains of course an open issue that 
deserves further study. 

The full 6D calculation we perform here is necessarily delicate because of the fact that 
the background geometry contains a conical singularity at the location of the 3-brane on 
which we intend to localize matter. This is the usual conical singularity associated with the 
tension of a codimension-2 defect. The orbifold symmetry imposed across the codimension-1 
brane creates a Z 2 reflection axis in the background conical solution which is further felt 
at the level of perturbations. Although there are many different ways of representing the 
background geometry, none is particularly straightforward for performing the perturbative 
analysis. To get a sense of the difficulty, consider the usual way of writing a conical defect, 



The angular variable runs from to 2ir, and 5 is the deficit angle proportional to the tension. 
The 3-brane is understood to be located at r = 0. The problem with this coordinate system 
is that the orbifold brane is located on two disjoint surfaces: = and = it. Even if the 
geometry is smoothed out at r = 0, it is clear that this metric will provide a poor description 
of physics on the codimension-1 brane near the codimension-2 brane. 
Alternatively, the geometry can be put in the form 





ds 2 = dr 2 + r 2 (l - 5) 2 d0 2 + rj^dx^dx" . 



(7) 



ds 2 = e^' 2) (dy 2 + d^ 2 ) + Vliu dx^dx u , with <f){y, z) = -5 ln{y 2 + z 2 ) . 



(8) 
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One advantage of this coordinate system is that the co dimension- 1 brane is now at z — 0. 
The drawback is clear, however: the conformal factor is singular at the location of the 
codimension-2 brane. One can deal with this by regulating the conformal factor, e.g., (f) — > 
—S ln(y 2 + z 2 + e 2 ), but this need to regulate the background geometry comes at the price of 
introducing extra terms in the perturbation equations which may or may not be negligible 
in the limit e — > 0. Although the final equations shall be independent of e, it is necessary 
to introduce the artifice of regularization to obtain them. 

Instead we will use a different form of the background metric, one that we have found is 
best suited for the analysis of perturbations: 

ds 2 = (1 + (5 2 (1 - e(y) 2 )) dz 2 + (dy + ^(y)dz) 2 + j^dx'W . (9) 

In this coordinate system, the codimension-1 brane is still at z — 0, and we work in the half- 
picture z > 0. The above metric is understood to be orbifolded at z — 0. The parameter 
(3 is related to the 3-brane tension. Meanwhile, e(y) is a regulating function with the 
following properties: e(oo) = 1, e{—y) = — e(y) and e'(y) = 28 e (y), where 5 e (y) is an explicit 
regularization of the Dirac delta function. The upshot of this coordinate system is that the 
induced metric on the codimension-1 brane is simply Minkowski space-time, and similarly 
for the induced geometry on the codimension-2 brane at z = y = 0. As we will see below, 
this greatly simplifies the analysis of perturbations. 

As always the presence of ghosts in a theory can be determined by checking the sign 
of the kinetic term of all propagating degrees of freedom in the action. In the case of 
gravitational theories, however, the inevitable presence of gauge modes and constraints 
complicates matters. Only after the gauge has been completely fixed and the constraints 
fully solved will the action reduce to the true physically propagating degrees of freedom 
whose kinetic terms can be inspected. As is well known, before the constraints have been 
solved the action may well contain wrong-sign kinetic terms, e.g. the famous conformal 
factor problem of Euclidean quantum gravity. 

One way to bypass this problem is to compute the coupling to a conserved source. Since 
gauge invariance enforces conservation of stress energy, any gauge degrees of freedom will 
not couple to a conserved stress-energy source, and, in particular, will not contribute to the 
single-graviton exchange amplitude. In the following we shall do precisely this to verify the 
absence of ghost degrees of freedom. 
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The paper is organized as follows. In Sec. II, we describe the background geometry (9), 
starting from the general ADM form of the 6D metric. In Sec. Ill, we study metric perturba- 
tions, focusing first on 4D vector and tensor modes. We then turn to 4D scalar perturbations 
in Sec. IV, and derive the lower bound on the tension. In Sec. V, we discuss couplings to 
external sources. We conclude with a summary of the results and discuss future avenues 
in Sec. VI. The Appendices include a review of the decoupling calculation (Sec. VII), an 
argument for why certain naively ghostly terms do not contribute to the imaginary part of 
the exchange amplitude (Sec. VIII), and an argument in (Sec. IX) for why one of the scalar 
modes, which at first sight looks problematic and ghostly, is in fact pure gauge and does 
not contribute to the exchange amplitude. Our notation throughout is as follows: A,B,... 
denote 6D space-time indices, a,b . . . are 5D indices, and u, . . . are 4D indices. The coordi- 
nates along the worldvolume of the 3-brane are denoted by x^ 1 , whereas the extra dimensions 
have coordinates y and z. The 4-brane is located at z — 0, and the 3-brane at y = z = 0. 



Our approach is to start with the full 6D action for the cascading setup written in the 
radial ADM formalism, and then perturb this action to quadratic order. We work in the 
"half (or reduced)-picture" , where we restrict ourselves to the z > side of the Z 2 -symmetric 
4-brane. Crucially this implies that all quantities evaluated at z = are understood in the 
following sense: (f>(z = 0) = lim 2 _> + <f>(z). 

In the ADM form, with z playing the role of a "time" variable, the 6D metric in the 
region z > is 



In what follows, covariant objects are constructed using the 5D metric g a b induced on 
constant-^ surfaces. The 6D part of the action, including the appropriate Gibbons-Hawking 
boundary term, is then 



II. CASCADING SETUP 



ds 2 = lAB dx A dx B = N 2 dz 2 + g ab (dx a + N a dz)(dx b + A 6 . 



■dz). 



(10) 




(11) 



with J dz = J Q °° dz, and where 




(12) 
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The full action is given by 

/+ f M A \ M 3 f 

+ f _ _ o d 4 x^ {^R, - \j , (13) 

where £ ma tter is the Lagrangian for the external sources which may be distributed in 6D, 
or localized on either brane. In what follows we will compute the coupling for an arbitrary 
conserved 6D, 5D or 4D source. 

A. Background Solution 

There are several ways to write the background solution. However, as discussed earlier, 
great care must be taken because of the singular nature of the geometry. Consider the 
following background metric 

ds 2 = ^ AB dx A dx B = (1 + (5 2 {l - e(y) 2 ))dz 2 + (dy + f3e(y)dz) 2 + rj^dx^dx" . (14) 

By comparing with (10), we can read off the background expressions for the ADM variables: 

N 2 = l + p 2 (l-e(y) 2 ); 
Ny = (3e(y); N» = 0; 

9ab = Vab- (15) 

The regulating function e(y) is such that e(oo) = 1, e(—y) = — e(y) and e'(y) = 25 e (y), 
where S t (y) is an explicit regularization of the delta function. The induced metrics on the 
codimension-1 (z = 0) and codimension-2 (z — y — 0) branes are both flat Minkowski space. 
The 6D Riemann tensor vanishes on the background for z > 0, as it should. For instance, 
since N 2 + N 2 = 1 + /3 2 , the Ricci scalar is given by 

Because this metric is orbifolded at z — 0, it follows from the Israel junction conditions 
that (14) correctly describes the metric of a codimension-2 object with tension localized on 
the orbifold plane. Explicitly, the Israel junction conditions are 

Mt{K5 a b - K\) =T%- MlG a b . (17) 

z=0 
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In the background, we have T a b = —X5(y) and G a b = 0. Meanwhile, in the gauge (10), the 
extrinsic curvature on the surface z = is 

K ab = ^ (d z g ab - V a N b - V b N a ) 

= ^ {d z g ab - g ac d b N c - g bc d a N c - N c d c g ab ) 
= -^V b . (18) 
At the background level, the Israel matching conditions (17) therefore read 

K 9 -A m - 5(y) (19) 

N v /l + ^ (1 _ e 2 (y)) M * W - U9j 

Since e'(y) = 25 e (y) — >■ 5(y), the extrinsic curvature does indeed encode the codimension-2 
source. The relation between the brane tension A and the required value for (3 requires some 
care, but can be determined unambiguously by integrating the extrinsic curvature, 

+V^. r dy , ^ m = [' W£) ^a.ctan,. (20, 



iv 7^ v /i + ^2 (1 _ e2(y)) y_ x v /i + /3 2 (1 _ e2(2/)) 

Thus the relation is 

A = 2Mg arctan /3 . (21) 

In particular, we recover the standard result that this static solution can only hold a maximal 
positive codimension-2 tension: A < ttM^ . (This is one-half of the usual condition because 
the Z 2 orbifolding projects out half of the space.) 



III. TENSOR AND VECTOR PERTURBATIONS 



To study perturbations, we perform a scalar-vector-tensor decomposition of the full 6D 
metric '-/ab — Iab + $1ab with respect to the 4D Lorentz group. The tensor and vector 
modes are straightforward as none of them couple to the positions of the branes and are given 
by a simple generalization of the massless scalar field. In particular, the tensor and vector 
sectors are manifestly ghost free, regardless of the value of the tension on the codimension-2 
brane. 
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Tensors: The metric for the tensor perturbations takes the form 

5<y A Bdx A dx B = hf^dafdx" , 
where rj^h^ = and d^h^ u = 0. The action for the tensor perturbations is given by 

S = ^ fdzj d 5 x (-Nd a h^d% u - i (£ n ^) 2 ) , (22) 

+ f / ^x(h^D 5 hl) + ^[ d*x(h^D^), 

° Jz=0 ° Jz=y=0 

where C n is proportional to the Lie derivative: C n X = d z X — N y d y X. Notice that C n and 
d y do not commute: 

d y C n X = C n dyX - dyNy dyX . (23) 

The tensor modes constitute 5 degrees of freedom, each of which satisfies the same equation 
which is manifestly ghost free since the action is only composed of a sum of positive kinetic 
terms. These five spin-2-degrees of freedom can be further decomposed into two helicity- 
2, two helicity-1 and one helicity-0 modes. The helicity-0 degree of freedom is expected to 
exhibit the Vainshtein mechanism as in usual massive gravity [6, 26, 27], DGP model [16, 28- 
31] and general theories of resonance graviton [3, 32]. 

Vectors: For the vector perturbations the metric takes the form 

<5 7 AB dx A dx B = 2(A^ + NyBfj) dx^dz + 2B^dy + (d^C u + d u C^ dx^dx" , (24) 

where the vectors A^, and are all transverse: d^A^ = d^B^ = d^C^ = 0. There is 
sufficient gauge freedom in the vector sector to set B^ = 0. The action for the remaining 
perturbations is then given by 

Af 4 C + ( 2 2 \ 

S = -± J d 6 x I {d y A,f + ^ (A, - CnC,) D 4 (A** - c n c») + 2Nd y c»n 4 d y cA 

+ _1 J d 5 x (2d y C»n,dyC,) + J d e xN (A,T» Z - C,d u Tn , (25) 

where T AB is an external 6D source, (which can include localized contributions on either 
of the branes). Note that neither vector field has a 4D source and so we expect that the 
vectors will be appropriately smooth at y — 0. With this in mind, we may neglect the mild 
^/-dependence of N, since its departure from N — 1 will be of zero measure in the integral. 
Then, on performing the following local field redefinition, 

= dyC^ ; A^ = A^ — C n C^ , (26) 
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the kinetic terms are diagonal and manifestly positive definite and so is the coupling to the 
source 

Skin = -f J d 6 x I^A^A^ + 2NC»n 4 cA + J tfx^CnkC^ (27) 
- J d 6 x NC^ (TM + N y T» z ) , (28) 



<7 

'-'sources 



where we used conservation of energy along the transverse direction V aT a ^ = 0. 

Notice that lives entirely in 6D, whereas has a kinetic term both in 6D and 5D. 
Crucially all the kinetic terms are positive, as required, confirming that there are no ghosts 
present in the vector sector. 



IV. SCALAR PERTURBATIONS 



The scalar perturbations are tricky because the positions of the brane transform as 
4D scalars. We choose to utilize the bulk gauge freedom to work in a gauge where the 
codimension-1 and -2 branes remain at fixed position, respectively z = and z = y = 0. 

At the perturbed level, the metric has seven 4D scalar modes: 

5 / yAB ( ix A dx B = h zz dz 2 + 2h zy dzdy + h yy dy 2 + 2h flz dx^dz + 2h /iy dx^dy 
+ {nritiv + d^dyw) dx^dx u . 

One can set one of these modes to zero by fixing the gauge while keeping the brane positions 
unperturbed. We use this freedom to set w = 0, such that h^ u = ni]^. This is an analogue 
of conformal Newtonian/longitudinal gauge used in cosmological perturbation theory. We 
are therefore left with six scalar degrees of freedom, which will be denoted by x? 0> V, a i T 
and Ti : 

h ab dx a dx b = irr] ab dx a dx b + Vdy 2 + 2<9 M r dx^dy 
h zz = 2N 2 <j ) + N a N b h ab + 2N a A b r ]ab 

h az = A b Vab + N b h ab , (29) 

where A a dx a = d a adx a + xdy- Note that A a is the perturbation in the shift vector, N a = 
N a + A a , while is the perturbation is the lapse function, N = N(l + <fi). 
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To quadratic order in perturbations, the full action is given by 



'quad 



' 6 / d 6 xC 6 / d 5 x7r(n4l / + 2n57r-2n 4 9j / r) 

J 4 J z=0 

d 4 a;^^7rn47r + A7r 2 ^ . (30) 

The kinetic term for 7r on the codimension-2 brane has manifestly the wrong sign. However, 
as we will see, the presence of brane tension A can change the sign of this kinetic term, 
after integrating out the non-dynamical degrees of freedom. Meanwhile, the mass term \n 2 
localized on the codimension-2 brane, which seemingly breaks the shift symmetry, will cancel 
out after integrating by parts a bulk contribution. 

To see how this works in detail, let us expand the 6D action to quadratic order, still 
focusing on the region z > 0. It is convenient to perform the following field redefinitions 

d N 4 

X = X ~ (C n + d y N y )r - 2^a + —C n d y n (31) 
n 4 a = U 4 a - 2C n 7r = -NSK^ (32) 
□ 4 f = n 4 r-2d y 7i. (33) 

In terms of these new field variables, the result for £ 6 at quadratic order is then 
£ 6 = ^X^x -n(v + 4tt- 2d y f + 2^<r) D 4 

-ld y N y C n {fU 4 f) + 3Nn(d 2 - D 4 )n - |(Or) 2 - 2%&d t * 2 

-in 4 cKAi - d y N y )(V + 4tt) - ^d y N y nD 4 a + ^D 4 ad y C n f + 3NnD 4 d y f 

3 - d N d N d 2 N 

--NttU 4 V + 2^^-aD 4 f - -j^-aTUa + d y NfD 4 (V + 4vr) . (34) 

The issue of the tt mass term in (30) can be immediately resolved. The last term in the 
second line of (34), when integrated by parts, generates a mass term for it on the 3-brane 
that precisely cancels that in (30): 



L 



-2M 6 4 / dz I dy^&d z n 2 = 2M 6 4 [ + dz [ dy (d z ^& ) n z - 2M 6 4 / dy 



N ./ ./ \ N 



dyNy 2 

N lz=o 



= 2Mt [ dy^n 2 = Avr 2 . (35) 

Jz=0 iV 

This cancelation makes manifest the shift symmetry of n. We now proceed to decipher the 
n kinetic term, which requires a few field redefinitions and integrating out non-dynamical 
fields. 
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A. Lagrange multipliers <p and V 

As it stands, (34) is linear in and V. Varying the action with respect to either of these 
therefore yields a constraint on the other degrees of freedom. We choose to vary with respect 
to and interpret the result as a constraint for V. Upon substituting V by its constrained 
value back into the action, we will see that the action picks up a quadratic term in a. 

The expression for V that follows from the variation is 

V = -4n + 2d y r- 2 -^a. (36) 

In the presence of an external source, there is an additional Tg 2 that enters in this equation; 
we shall take this into account in Sec. V when we compute the relevant couplings — see 
Eq. (49). Substituting this expression for V in the 6D action, we can simplify the last two 
lines in (34) as follows 

A = -in 4 cr(£ n - d y N y )(V + 4tt) - -^d y N y nU 4 d + -^U 4 dd y C n f + ?>NnU 4 d y f 

--NttD 4 V + 2^Si^aD 4 f - ^£aD 4 a + d v NfD 4 (V + 4tt) 

1 - / 2d N \ 3 2 

= ~^n 4 a(C n - d y Ny) [2d y f - ~r~°} - ^d y N y nn 4 a + ^D 4 ad y C n f 

+3Nnn 4 d y f - ^A%D 4 (-4vr + 2d y r - ?^r) + 2^^W 
<9 2 iV _ / d N \ 

-^r<W + d y Nfn A [2d y f - 

2 - (d N \ d 2 N 

= QNnD^ - d 2 y NfU 4 f + J 4 a(£ n - d y N y ) I -j^a) - -^-aU^a , (37) 

where in the last step we have used (C n — d y N y )d y f = d y (C n f) from (23). We can further 
simplify this expression by integrating by parts: 

r+ r+ ( _ 9 (d N \ d 2 N 

/ d 6 x A = d 6 xi 6NttD 4 tt - d 2 y NfU 4 f + -^U 4 a{C n - d y N y ) f -^-a J - ~^ 4 a 

/ + - f f) N 

d 6 x (6A7rD 4 7r - ^7VfD 4 f) - / d 5 x^/aD 4 a , (38) 
J z =o N 



where the last follows because N 2 +N 2 = 1+/3 2 = constant. Substituting (38) back into (34), 
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we note that the quadratic terms in r simplify through integration by parts: 

J 2 ' 

2 ./ y \ N ~ 



dz[%N+^£ n ) (fD 4 f) = -\ [ + dzd y 



1 TU4T 



9 "%n t f 



N 



J 

Jz=0 



d N 
N 





(39) 



Combining all of these results, the 6D part of the quadratic action reduces to 
f + d 6 xC 6 = J + dzd 5 x (JfiXEkX + 3A%D 5 7T - |(Or) 2 

Similarly, we must also substitute (36) for V into the 5D part of (30): 

3 ^ 3d v N v ^ 



£-5 — 2 #5-^5 



(40) 



(41) 



quadratic 

As we will see, the above kinetic mixing term between 7r and a is what provides an additional 
4D kinetic term for tt. 

Putting everything together, and using the background Israel Matching condition (19), 
M^dyNy/N = X5(y), the complete action for the scalar modes is given by 
M K 4 



S, 



quad 



d b x 



±xn,X + z(Nnn 5 n-l(£ n n) 2 



+ Ml I d 5 x^7rD 5 7r 



+ Mi 



2=0 



/ 

Jz=y=0 



3 
2' 



d A x 



3^ A/ 1 ^ 3 ^ 

4 2M| \ N 2 m 6 N 



• (42) 



B. a degree of freedom 



The last step consists of diagonalizing the kinetic matrix for ir and a. This is achieved 
by the following change of variable 

„ _ 2N 



a = a 



7T . 



(43) 



in terms of which the action becomes 



'S'quad 


Ml 
2 




+ 


Ml 


[ d 5 

lz=0 


+ 


Mj 


/ 

l z =y=0 



ixD4X + 3^iV7rn 5 7r-i(£ n 7r) 2 



d 4 x 



3A 



4 \2m 2 M 2 



1 7iD 4 n + 



A 



2 , rD 4 f - ^*n 4 * 



1 

N 2 ' 



• (44) 
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This is our main result. Through a series of field redefinitions, and after integrating out 
auxiliary fields, the resulting kinetic term for it on the codimension-2 brane acquires a A- 
dependent contribution: 

^^(dpr 1 )^- (45) 

Thus the ir mode is healthy as long as the tension is larger than the lower bound 

A mi „ = 2 -Mlm\. (46) 



V. COUPLING TO AN EXTERNAL SOURCE 

It is immediately apparent from (44) that n, x and f are manifestly unitary. However, 
a appears as a ghost and would seem to present a serious concern. (This issue never arose 
in [9] since this mode disappears in the decoupling limit.) As already discussed, though, it is 
common in the analysis of gravitational systems to find apparently ghostly degrees of freedom 
that nevertheless decouple when a physical, gauge-invariant calculation is performed. In this 
case we will see that neither a nor r couple to conserved matter, and consequently these 
modes are pure gauge, at least to this order in perturbation theory. To see this, let us 
introduce the matter content in the action 

dz J d 5 xC mattcI = fdzj d 5 x^h AB T AB : (47) 

where we consider a 6D stress-energy tensor in all generality. In particular the total stress- 
energy T^ B can include terms localized on each of the branes: 

T 6 AB = T AB + N-H+iz) % ab 5X + N-'5(y)5 + (z) 7T W > ( 48 ) 

where 5+(z) is the delta function normalized on the half-line: J °° dz 5+(z) = 1. However, 
it is not necessary for us to distinguish between these different contributions since, in any 
case, any source must be covariantly conserved in a 6-dimensional sense. In Appendix B 
we discuss a subtle issue having to do with conservation of 5-dimensional external sources 
of the above form that appears to arise because of the singular nature of the background 
geometry. 

In the presence of an external source, the constraint equation obtained by varying with 
respect to gets modified from (36) to 

2d N N 2 

V^-4. + 2^--^+_T«. (49) 
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On substituting this back into the action, everything that is linear in drops out, including 
part of the matter action. The remaining matter contribution which sources (44) is 

d v N v 



a 



matter 



= -No 

-Nf 
3N 2 



BaT Az + " y ^' y (T : 
6 N 2 { 6 



mi 



NyTD 

NA,N 



d A T 6 Ay + N y d A T Az - i^p(rf + N yT p 



-IT 



2m§ 



d A Tt 



N 2 
ON 



+ Nx Tf + N y T* 



+ Nn[2T^-^ AB T AB 



3^ 
2 



N d N 

a a i Ft Ft T m V V V a 



N 2 



N 3 



1 



Ml 6 4 + ie 



[Tyy + 2N V T»* + N 2 T* Z ) 



-N&vfT Az 



Nf 



v A J 6 ly y w A 1 6 



+ Nx 



3N 



nVfT Az + 2N—d^ { H ) T 6 A » + Ntt 

6 ^4 



T yz 4- N T z: 
3 



2T* - ~1 AB T AB - 2^ 



rpfMU 



N 3 

rpz; 

Ml 6 



1 



(Tf + 27V,Tf + N 2 y Tl z ) , 



(50) 



where Tf M 



□ 4 + ?e 

r^Tf , and is the 6D covariant derivative. In performing several integra- 
tion by parts we have made the assumption that T 6 Zj4 | 2=0 = 0, which ensures that there is no 
external force on the orbifold plane. This condition is only necessary in the 'half-picture' we 
are using here. In the doubled picture these terms are canceled by identical contributions 
of opposite sign on the other side of the orbifold plane. 

We see that both modes a and f couple only to the transverse part of the stress-energy 
tensor. For conserved matter, V^Tf 8 = 0, both a and f decouple, and the matter sector 
only sources x an d 7r 



■^matter 

N 3 



Tf + N y T t 
1 



+ Ntt 



-iAB 



rpzz 

Mr 6 



2T 6 ; - ~1 AB T A 
(Tf + 2N y Tl z + N 2 T ZZ ) . 



c^ d^dy rplXV 

1 T37 6 . 



(51) 



i 6 n 4 + ie 

The sign of the kinetic term of a in the final action (44) is therefore irrelevant to this analysis, 
since this mode is clearly pure gauge and brings no instability. Furthermore, we notice that 
X only couples to 6D matter in the bulk, and thus n is the only scalar mode to be excited 
when including an external source on the branes only. Since x is n °t a true 6D degree of 
freedom, we may as well integrate it out, leaving only an action for n with additional source 
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contributions: 



+ f d e x Nn [27ft - \^ AB Tt B - 2^T 6 ^ 

+ / + d " x ( ^nTv^ " ^nTv^) " (52) 
y m 6 V U 4 -Me U 4 + ze / 

Now the additional source contributions may appear to have a 4-dimensional pole, but this 
is in fact not the case. Such terms arise whenever 6D perturbations are decomposed into 4D 
scalar- vector-tensors, and appear non-local because of the inherent non- locality of the split. 
In Appendix A we demonstrate that these do not contribute a pole to the single-graviton 
exchange amplitude. 

When considering matter on the 3-brane only, T^ B = N~ 1 5(y)5 + (z) TfS^S^ , with 
d^T^ v = 0, the matter Lagrangian is simply 

^malter = ■ (53) 

This agrees with the flat space-time limit, as well as with the decoupling limit. Meanwhile, 
the Lagrangian for 5D matter is 



£ (5) = -it 

■'-matter ^ 



77W4^7T , (51) 
u 4 



where we used d ii d v T§ v = dyl^ v , which follows from 5D stress-energy conservation. 



VI. DISCUSSION 



In this paper we have demonstrated the perturbative unitarity and stability of the 6D 
cascading gravity framework. This puts the claim that this model is a consistent infrared 
modification of gravity on a firm footing. While we cannot rule out from the previous 
analysis the existence of an non-perturbative instability, we at least expect a finite range 
of energies for which this model can be understood as consistent effective field theory. A 
potentially profound result of this analysis is the confirmation of the result, first obtained 
in the decoupling limit [9] , that unitarity requires a minimum tension for the codimension-2 
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brane localized on an orbifold plane, as soon as we introduce a localized kinetic term for 
the graviton, i.e. a 4D induced gravity term. It will be interesting to explore whether a 
similar bound on the tension applies to higher-codimension branes, and whether there is any 
fundamental reason for this bound, or whether it is possible to relax this condition by other 
means (see [10, 24]). In any case, our results establish this model as a consistent framework 
within which to explore the phenomenology of infrared theories of modified gravity and to 
confront its predictions against cosmological observations. 
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VII. APPENDIX A: THE DECOUPLING ARGUMENT 



We briefly review the analysis of [9] using the decoupling limit. (Note the slight notational 
differences compared to [9]. In particular, we have interchanged the roles of h ab and h ab 
below.) In the decoupling limit M 5 ,M 6 — > oo with A 6 = (mgMf) 1 / 7 fixed, the action (1) 
reduces to a local field theory on the orbifold brane, describing 5D weak-field gravity and a 
self- interacting scalar field it: 



oJordan 
'-'decoup 



-V 



+ 



d 5 x 
d 4 x 



l -h a \£h) ab + Tr V ab (£h) ab - ^L(d a 7r) 2 n 5 7r 



Mi 



h^{£h), v ■ - 2 h" • 



(55) 



where (£h) ab = -\ (D 5 h ab - i] ab D 5 h - d a d c h cb - d b d c h ac + r] ab d c d d h cd + d a d b h) is the lin- 
earized Einstein tensor in 5D, and (Eh)^ that in 4D. All other interactions in (1) are 
suppressed by powers of 1/M 5 , 1/M 6 and therefore drop out in the decoupling limit. 
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With the action written in terms of the physical metric h ab , as in (55), we notice that n 
is non- minimally coupled to the linearized Ricci scalar, i.e., h ab is the Jordan-frame metric. 
The kinetic matrix for the metric and 7r can be diagonalized as usual by shifting to the 
Einstein frame variable, h ab = h ab — Ttr] ab . Expressed in terms of h ab , the decoupling action 
becomes 



oEinstcin 
'-'decoup 



d 5 x 



r / (I' 

Jy=0 



1 + 



9 



16m; 



■□ 5 7T 



M 2 1 ~ 1 

.-±h^(£h)»„ + -h^ + -nT 



(56) 



As usual, the transformation to Einstein frame generates a conformal coupling of 7r to matter. 

In the presence of tension on the 3-brane, T^ u = —\r]^ u , the solution to the equations of 
motion is 



7T 



(0) 



A 



3M| 



\y\ ; 



h\2 



A 



3M| 



\y\vvw ; 



u yy 



0. 



(57) 



These combine to give a flat Jordan-frame metric, h 



(o) 



h y ^l + TT^r]^ = 0, consistent 



with the fact in 6D that the codimension-2 source generates a deficit angle but leaves the 
geometry Riemann flat. 

To disentangle the different degrees of freedom, let us define 



h ab dx a dx b = h^dx^dx" + Vdy 2 + 2V lx dx»dy . 



(58) 



Substituting into the action (56) gives 



oEinstein 
'-'decoup 



d b x 



Ml 



2V(d»d% u - U 4 h {A) ) + W»{d y dJi {A) - d y d% u ) 



+h^n 5 v - hpsh) + 2^(4)^9" V + 2(^ v) 2 - f^) 



3M| 



(9a 



7T 



+ 



/ 



16ml 



□ 5 7T 



d 4 x 



M 2 1 ~ 1 

-fh^{Sh), v + -h^T, u + -nT 



(59) 



with F M y = d^V v — d v V^. Concentrating on scalar perturbations only, we may choose a gauge 
in which the Einstein metric is pure trace: h^ u = ^h^rj^ v . However, now it is clear that V 
and the scalar part of are Lagrange multipliers which enforce the constraint O^h^ = 0. 
Thus in the scalar sector we have O^h^u = O^r)^ and O^h^ = 0. Note that these equations 
are consistent with the above background solution, at least for \y\ > 0. The action therefore 
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reduces to the simple form (neglecting the vector part of which does not couple to 4D 
matter) 

9 



on 
^decoup 



d*X 



1 + 



16ml 



□ 5 7T 



+ 



d*X 



3M 2 ^ 1 m 
— - — ttU^ty H — ixT 
4 2 



(60) 



Expanding (60) to quadratic order in perturbations around the background (57), we find 
that the i\ cubic term yields a contribution to the kinetic term of the perturbations that is 
localized on the 3-brane: 

/9A 
d 5 x^(d,n) 2 5(y), (61) 

where fc now denotes a perturbation around the background profile tt^°\ Combining with 
the rest of (56), the quadratic action is thus given by 



C7T 

^dccoup 



/ 



d 5 x 



7T~ \ d a n ) 



+ 



d x 



3M| 



3A 



- 1 7rn 4 7r + -n5T 



2m 2 R M? ~ I ' 2 



(62) 



This is now manifestly equivalent to the decoupling limit of our final result, reproducing the 
same minimum bound on the tension. 



VIII. APPENDIX B: NO EXTRA POLE IN TT AMPLITUDE 

In the final action for 7r, given in (52), we obtained the following term 

W = J +dz J d5 *5 ( T " D7T¥ r " " rf nrT¥ r « ') ■ < 63 > 

Superficially this looks like a contribution to the TT amplitude which has a 4 dimensional 
pole and will in general be ghostly. However, since this term already arises in Minkowski 
space-time it is clear that this cannot be the case. This non-local term arises because the 
scalar-vector-tensor decomposition is inherently nonlocal. 
Absence of ghosts is equivalent to the condition that 

lm[W] > . (64) 

Since the imaginary part comes entirely from the pole at D 4 = 0, it is sufficient to consider the 
form of the stress-energy on the pole surface, i.e., in 4-momentum space with k 2 = 0. Since 
we are looking at scalar perturbations it then follows that on this surface d^T^ = d^T^ = 0. 
Now conservation of energy implies 

d z (NT^) + d y {NTD + 9 -^Tl V = ; (65) 
d z (N 2 T* y ) + d y (N 2 T yy ) = . (66) 
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We can solve these equations as follows 

Tyy = N-*d 2 z U, Tt = -N-*d z d y U, T** = ld y {N-%U)-^d z U. (67) 

Computing the relevant piece, and using the fact that the Z 2 symmetry implies T^ = 0, 
then d z U = at z = 0, and so we find when expressed in momentum space 

J + dzdyN* (T^(-k)Tyy(k) - Tt(-k)Tt(k)) 

= j + dz&y (-d, (N^dyU (-k)d z dyll (k)) - ^d z \d z U{k)\^j 

= (N-%U(-k)d z d y U(k)) z=0 + ^\d z U(k)\ 2 z=0 = 0, (68) 

where k is the on-shell 4-momentum. In other words, conservation of energy and the orbifold 
condition ensures that there is no anomalous 4D pole in the propagator. 



IX. APPENDIX C: NO POLE FROM a IN TT AMPLITUDE 

Although the stress-energy on the 4-brane is automatically conserved in 5D and 6D, it is 
worth stressing that the presence of the singularity at y = z = imposes an extra constraint 
for the allowed matter on the orbifold brane. Indeed, the 6D conservation of energy condition 
applied on the 5D stress-energy implies 



^% Az S + (z)] = T^ B % AB N-H + {z) 



= ^V y S + (z) = ^LpT?>6(y)6 + (z) = . (69) 

In the background gauge (10), the 6D conservation of energy imposes the ^-component 
of the 5D stress-energy to vanish at the singularity T^ y = 0. If this is not so, then 

y=z=0 

this must be compensated by a nonzero contribution to T zA . This is a surprising condition, 
and only really becomes of importance because of the need to regulate the background 
geometry. Although at first sight it would seem to forbid putting a very small tension on 
the codimension-1 brane, it is very likely that in practice in the presence of a codimension-1 
tension the background will be sufficiently modified that this condition is removed. 

If we do not make the assumption that T yy = at y = z = 0, then on coupling to the 
metric induced on the orbifold brane, it would appear that a, which decouples from 6D 
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conserved matter, is nevertheless sourced by a term of the form 

/ ^ - 5Aj7?» = - / d%-^7?» . (70) 

Since a enters with the wrong sign kinetic term, such a coupling would seem to indicate a 
violation of unitarity. Integrating out a this piece generates a term of the form 

/ ^A^u^' - 1 ^ar^' < 71 > 

However, following the same reasoning as in Appendix B, since the unitarity violation comes 
entirely from the pole, we may look at the implications of the 5D conservation equations on 
the pole surface k 2 = 0. Since we are looking at scalar perturbations it again follows that 
9^7^^ = on this surface, and so we have 

dylT = • (72) 
For a localized source we may integrate (71) by parts to rewrite it as 

Im (| ^y)\V y ^—V y ) = "Im (/ d ^<y)V V ^-dyV y ^ = . (73) 

Thus there is no pole contribution to the TT amplitude coming from <j, further confirming 
that its apparent wrong sign kinetic term does not indicate the presence of a ghost. 
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